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Movements  of  upper  trophic  level  predators  in  an  open  ocean  environment  should 
be  driven,  in  part,  by  the  distribution,  density,  and  movements  of  their  prey.  Sur¬ 
veys  have  shown  that  cetacean  densities  are  higher  closer  to  shore  around  the  main 
Hawaiian  Islands  than  in  offshore  waters  (Barlow  2006),  presumably  reflecting  in¬ 
creased  productivity  or  spatial  and  temporal  predictability  of  prey  associated  with 
island  effects  (Baird  et  al.  2008 a).  A  number  of  high  trophic  level  pelagic  species  have 
been  shown  to  concentrate  around  and/or  use  mesoscale  eddies  as  foraging  habitat 
(e.g. ,  Davis  et  al.  2002,  Seki  et  al.  2002,  Bakun  2006,  Polovina  et  al.  2006,  Yen  et  al. 
2006).  The  islands,  and  their  interaction  with  winds  and  currents,  create  a  complex 
system  of  eddies  that  may  also  concentrate  some  prey  types  farther  offshore  (Seki 
et  al.  2002),  but  whether  island-associated  cetacean  populations  use  these  offshore 
eddy  systems  for  foraging  habitat  is  unknown. 

Melon-headed  whales  (. Peponocephala  electro)  are  a  tropical  delphinid  that  primarily 
inhabits  open-ocean  waters  and  is  known  to  feed  on  mesopelagic  squid  and  fish 
(Best  and  Shaughnessy  1981,  Clarke  and  Young  1998).  In  Hawaiian  waters,  these 
whales  are  encountered  relatively  infrequently.  Around  the  main  Hawaiian  Islands, 
two  populations  have  been  identified:  a  “main  Hawaiian  Islands”  population  and  a 
population  resident  to  the  island  of  Hawaih  (Aschettino  et  al.  2011).  As  part  of  a 
multispecies  study  examining  movements  and  habitat  use  of  cetaceans  around  the 
main  Hawaiian  Islands  (e.g.,  Schorr  et  al.  2009,  Baird  et  al.  2010),  individual  melon¬ 
headed  whales  from  the  “main  Hawaiian  Islands”  population  have  been  instrumented 
with  satellite  tags  for  various  purposes  (e.g. ,  Baird  et  al.  2008 b),  including  to  assess 
their  movements  in  relation  to  mesoscale  eddies,  to  determine  whether  they  use 
offshore  eddy  systems  as  foraging  habitat. 

To  test  the  hypothesis  that  movements  of  melon-headed  whales  leaving  the  vicin¬ 
ity  of  the  main  Hawaiian  Islands  will  be  influenced  by  offshore  eddy  systems,  we 
analyzed  satellite  telemetry  data  from  10  melon-headed  whales  in  conjunction  with 
output  from  the  Hybrid  Coordinate  Ocean  Model  (HYCOM)  Hawaih  regional  model 
(http :  /  /  apdrc .  soest  .hawaii .  edu/REGMOD) . 

Details  on  the  satellite  tags  and  deployment  system  were  presented  by  Schorr  et  al. 
(2009)  and  are  only  briefly  discussed  here.  The  Low  Impact  Minimally  Percutaneous 
External-electronics  Transmitter  (LIMPET)  tags  (Andrews  et  al.  2008,  Schorr  et  al. 
2009,  Baird  et  al.  2010)  deployed  were  location-only  SPOT5  (Wildlife  Comput¬ 
ers,  Redmond,  WA)  Argos-linked  platform  transmitter  terminals,  which  attached 
to  the  dorsal  fin  using  two  4.5  cm  long  medical-grade  titanium  darts  with  back¬ 
wards  facing  barbs  that  acted  as  anchors  for  the  darts  (see  Andrews  et  al.  2008). 
Tags  were  deployed  remotely  using  a  Dan-Inject  JM  Special  25  (Bprkop,  Denmark) 
pneumatic  projector  with  a  modified  arrow  to  hold  the  tag  in  flight  at  a  range  of 
3-6  m. 

Age  class  of  tagged  individuals  was  characterized  in  the  field  based  on  relative  body 
size  as  “adult”  or  “juvenile.”  Tagged  individuals  were  photographed  and  photos  were 
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compared  to  a  photo-identification  catalog  (see  Aschettino  et  al.  201 1)  to  determine 
population  identity. 

We  used  the  Douglas  Argos-Filter  (available  at  http://alaska.usgs.gov/science/ 
biology/spatial/douglas.html),  version  7.06,  to  assess  locations  for  plausibility  us¬ 
ing  two  independent  methods;  distance  between  consecutive  locations,  and  rates 
and  bearings  among  consecutive  movement  vectors  (see  Schorr  et  al.  2009).  The 
user-defined  variable  of  maximum- redundant  distance  (temporally  near-consecutive 
points  within  a  defined  distance  that  are  kept  by  the  filter)  was  set  at  3  km,  maximum 
rate  of  movement  was  set  to  20  km/h,  Argos  class  2  and  3  locations  were  always 
retained,  and  a  rate  coefficient  of  25  was  selected.  The  rate  coefficient  algorithm  takes 
into  account  that  the  farther  an  animal  moves  between  locations,  the  less  likely  it 
is  to  return  to,  or  near  to,  the  original  location  without  any  intervening  positions, 
creating  an  acute  angle  characteristic  of  typical  Argos  error.  Results  discussed  below 
use  only  locations  that  passed  these  filtering  parameters.  From  these  locations,  only 
the  position(s)  with  the  best  location  class  for  each  day  were  used.  On  days  when  there 
were  multiple  positions  with  the  best  location  class,  the  positions  were  averaged  to 
obtain  one  daily  position. 

Tracks  were  broken  into  nearshore  and  offshore  components.  Offshore  components 
began  on  the  first  day  when  all  subsequent  positions  were  in  waters  deeper  than 
3,000  m.  A  depth  of  3,000  m  was  used  based  on  the  steeply  sloping  bathymetry 
around  the  Hawaiian  Islands  and  because  it  resulted  in  the  clearest  separation  be¬ 
tween  track  segments  that  remained  offshore  and  segments  that  returned  to  shore. 
Bathymetry  data  are  described  in  Smith  and  Sandwell  (1997). 

Daily  sea  surface  height  (SSH)  and  currents  with  0.08°  resolution  from  the 
HYCOM  Hawaih  regional  model  were  examined.  Horizontal  gradients  in  both 
variables  were  used  to  identify  eddies.  In  this  study,  HYCOM  data  were  used 
in  place  of  AVISO  satellite  data  as  the  former  are  available  daily  and  the  latter 
only  as  weekly  means  (http://  www.aviso.oceanobs.com/en/data/products/sea-surface- 
height-products/global/index.html).  Additionally,  AVISO  SSH  and  current  data  have 
only  a  0.3°  resolution.  The  finer  spatial  and  temporal  resolution  offered  by  the 
HYCOM  output  allowed  for  more  detailed  analyses  of  the  melon-headed  whale- 
selected  habitat.  HYCOM  assimilates  satellite  altimetry  and  therefore  should  repre¬ 
sent  the  mesoscale  features  observed  in  the  AVISO  SSH  data  well.  To  confirm  this, 
the  daily  HYCOM  SSH  data  were  regridded  to  AVISO’S  0.3°  resolution  and  weekly 
means  calculated  for  the  time  and  area  spanned  by  the  recorded  melon-headed  whale 
tracks.  These  regridded  means  were  well  correlated  with  the  AVISO  data  (r  =  0.88), 
as  were  the  regridded  daily  HYCOM  values  with  AVISO  weekly  means  (r  =  0.85). 

A  Kolmogorov-Smirnov  (KS)  test  was  used  to  investigate  whether  melon-headed 
whales  selected  a  significantly  different  environment  than  the  surrounding  area. 
As  discussed  by  Kobayashi  et  al.  (2008),  the  KS  test  compares  two  cumulative 
distributions  by  using  the  maximum  vertical  distance  between  the  two  distributions 
as  the  test  statistic.  If  the  distribution  of  the  specific  oceanographic  variable  along 
the  track  is  significantly  different  than  that  of  the  surrounding  area  (defined  below), 
it  can  be  inferred  that  the  animal  is  selecting  its  habitat  based  on  the  variable  being 
examined. 
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KS  tests  were  performed  with  both  horizontal  surface  current  magnitude  and 
surface  vertical  velocity.  Horizontal  current  magnitude,  |V|,  is  defined  as: 

|vi  =  Vl«2  +  u2i  (i) 

where  u  and  v  are  the  zonal  and  meridional  components  of  the  surface  current, 
respectively.  Distributions  of  these  variables  coinciding  with  the  offshore  track  posi¬ 
tions  were  compared  to  distributions  of  those  in  the  surrounding  area  within  defined 
bounding  limits.  Each  daily  track  position  was  compared  with  the  HYCOM  output 
for  the  corresponding  date,  and  these  daily  data  were  aggregated  to  form  the  track 
and  surrounding  area  distributions  used  in  the  KS  tests.  The  track  distributions 
include  only  data  from  grid  cells  encompassing  daily  track  positions.  Surrounding 
area  distributions  include  all  remaining  data  within  the  bounded  surrounding  area. 
The  latitudinal  limits  of  the  bounded  surrounding  area  were  chosen  such  that  they 
included  the  area  where  the  melon-headed  whales  could  have  traveled  and  were  set 
as  follows: 


MaxLatsA—  MaxLatTrack+0.5  (MaxSpan)  (2  a) 


MinLatsA=  MinLatTrack  —  0.5  (MaxSpan),  (2b) 

where  MaxSpan  represents  the  greater  of  either  the  total  latitudinal  or  longitudi¬ 
nal  distance  spanned  by  the  offshore  portion  of  the  track.  The  same  method  was 
used  to  determine  the  minimum  and  maximum  longitude  of  the  surrounding  area. 
Longitudinal  distances  were  calculated  in  kilometers  at  the  mean  latitude  for  the 
offshore  portion  of  the  track.  To  test  the  sensitivity  of  the  KS  test,  values  of  both 
0.2 5 (MaxSpan)  and  0.7 5 (MaxSpan)  were  tested  as  well. 

Tracks  ranged  from  7  to  24  d  (x  =  1 5  d;  SD  =  5  d)  and  from  172  to  1 ,140  km  (x  = 
627  km;  SD  =  319  km)  in  length.  All  tagged  whales  were  classified  as  adults.  Five 
of  the  whales’  tracks  (tagged  in  four  different  groups)  contained  offshore  components 
(Table  1,  Fig.  1),  while  the  remaining  five  were  entirely  nearshore.  The  offshore  track 
components  ranged  from  4  to  16  d  (x  =  9d;  SD  =  5  d)  and  from  160  to  819  km 
(x  =  463  km;  SD  =  283  km)  in  length.  Although  two  of  the  five  whales  were  tagged 
within  the  same  group,  they  remained  together  only  when  nearshore  and  separated 
for  the  offshore  portion  of  their  tracks  (Tracks  1  and  2,  Fig.  1). 

Two  tracks  near  the  edges  of  cyclonic  eddies  were  characterized  by  greater  |V| 
than  the  surrounding  area  (Tracks  4  and  5,  Fig.  2A— D,  Table  1),  while  the  opposite 
was  true  of  the  track  near  the  center  of  an  anticyclonic  eddy  (Track  3,  Fig.  2E-F, 
Table  1).  Sensitivity  tests  using  both  0.25(MaxSpan)  and  0.75(MaxSpan)  also  yielded 
significant  results  (1.6  x  10-3  <  P  <  4.7  x  10-2)  for  tracks  3,  4,  and  5  (0.47  <  P  < 
0.9)  for  tracks  1  and  2).  Thus,  we  feel  bounding  the  KS  test  area  by  0.5 (MaxSpan) 
adequately  represents  the  area  the  melon-headed  whales  could  have  selected  for 
transit.  KS  tests  on  vertical  currents  were  insignificant. 
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Table  1.  Details  of  all  melon-headed  whale  tracks  with  an  offshore  segment.  Numbers  in 
parentheses  indicate  start  dates  and  lengths  for  the  offshore  segment  of  the  track.  The  mean 
horizontal  current  magnitude  |  V  |  is  shown  for  the  offshore  track  and  bounded  surrounding 
area,  excluding  the  offshore  track,  used  in  the  KS  test. 


Mean  |V|  (m/s) 


Track 

Start  date 
(offshore) 

End  date 

Length  (km) 
(offshore) 

%  Time 
offshore 

Track 

Excluding 

track 

KS  P -value 

1 

19  Apr  08 
(02  May  08) 

12  May  08 

1,140  (551) 

46 

0.26 

0.23 

0.66 

2 

19  Apr  08 
(23  Apr  08) 

28  Apr  08 

557  (187) 

40 

0.18 

0.23 

0.59 

3 

25  Jun  08 
(29  Jun  08) 

14  Jul  08 

889  (819) 

80 

0.16 

0.23 

1.8  x  10“2 

4 

02  Jul  08 
(12  Jul  08) 

20  Jul  08 

1,061  (599) 

47 

0.29 

0.17 

3.2  x  10“2 

5 

10  Dec  08 
(20  Dec  08) 

24  Dec  08 

543  (160) 

33 

0.36 

0.22 

5.3  x  10-3 

Both  statistical  and  visual  examination  of  the  offshore  portion  of  three  melon¬ 
headed  whale  tracks  show  that  melon-headed  whales  around  the  main  Hawaiian 
Islands  that  move  offshore  occupy  the  edges  of  cold  core  cyclonic  eddies  and  the 
centers  of  warm  core  anticyclonic  eddies,  presumably  for  foraging  purposes  (Fig.  2). 
Visual  inspection  confirmed  that  the  melon-headed  whale  tracks  were  indeed 
coincident  with  these  specific  eddy  regions,  while  statistical  analysis  was  used  to  test 
whether  the  currents  in  these  regions  were  significantly  different  than  those  in  the 
surrounding  area.  Convergence  zones  in  cyclonic  and  anticyclonic  eddies  are  marked 
by  strong  and  weak  currents,  respectively,  eliminating  the  possibility  that  these 
regions  were  instead  chosen  for  physical  characteristics  such  as  a  preferred  surface 


Figure  1 .  Tracks  of  all  10  melon-headed  whales  from  the  main  Hawaiian  Islands  population 
with  nearshore  portions  in  gray  and  offshore  portions  in  black.  Downward-pointing  triangles 
represent  tag-deployment  positions  while  upward-pointing  triangles  represent  final  recorded 
track  positions.  The  five  offshore  tracks  are  numbered  near  their  final  recorded  track  positions. 
The  rectangle  represents  the  boundaries  of  the  Ni'ihau  Box. 
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Figure  2.  Melon-headed  whale  tracks  4  (A),  5  (C),  3  (E),  and  1  (G)  and  HYCOM  SSH 
(color)  and  surface  currents  (vectors)  for  13  July  2008,  21  December  2008,  8  July  2008,  and 
9  May  2008,  respectively.  The  nearshore  segment  of  the  track  is  represented  by  a  dotted  line, 
offshore  by  solid.  The  cyclonic  eddies  with  which  the  tracks  4,  5,  and  1  are  associated  are 
indicated  by  low  SSH  (<65  cm)  and  cyclonic  rotation  in  the  surface  currents.  The  anticyclonic 
eddy  with  which  track  3  is  associated  is  indicated  by  high  SSH  (>75  cm).  The  red  circle 
designates  the  position  of  the  melon-headed  whale  on  the  day  coinciding  with  the  HYCOM 
output.  Triangles  as  in  Figure  1 .  Area  shown  in  (A),  (E),  and  (G)  represents  area  enclosed  by  KS 
test  bounding  limits.  Rectangle  in  (C)  represents  area  enclosed  by  KS  test  bounding  limits. 
Box  in  (G)  indicates  the  Ni'ihau  Box.  Both  track  (black  solid)  and  nontrack  surrounding  area 
(blue  dashed)  cumulative  distributions  used  in  the  KS  test  are  shown  for  tracks  4  (B),  3  (D), 
3  (F),  and  1  (H). 
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current.  Both  warm  and  cold  core  eddies  are  continually  recurring  features  in  the 
lee  of  the  Hawaiian  Islands.  Wind  stress  from  easterly  trade  winds  is  intensified  by 
the  islands’  topography,  leading  to  the  formation  of  eddies.  These  eddies  then  propa¬ 
gate  northwestward  (Patzert  and  Wyrtki  1974,  Vaillancourt  etal.  2003).  To  the  south 
of  the  islands,  the  westward-flowing  north  equatorial  current  generates  and  sustains 
westward  and  northwestward  moving  eddies  as  well  (Holland  and  Mitchum  2001). 
The  melon-headed  whale  tracks  also  indicate  a  potentially  favorable  habitat  south¬ 
west  of  the  Hawaiian  island  of  Nihhau  (Fig.  1),  although  the  factors  contributing  to 
this  area’s  desirability  are  unclear. 

Two  melon-headed  whale  tracks  were  characterized  by  their  association  with  the 
edges  of  cyclonic  eddies;  one  spending  4  of  its  5  offshore  days  in  this  edge  region 
and  another  6  of  its  9  offshore  days  (Tracks  4  and  5,  Fig.  2 A— D).  Visual  inspection 
confirms  that  both  tracks  associated  with  cyclonic  eddies  were  associated  with  the 
outer  edge  region.  This  edge  region  is  marked  by  greater  currents  than  the  eddy 
center,  as  can  be  seen  in  Figure  2.  Additionally,  this  edge  region  is  characterized  by 
convergence  of  nutrients  and  phytoplankton  upwelled  in  the  eddy’s  divergent  center. 
Several  studies  identify  this  convergence  zone  as  a  fruitful  foraging  ground  for  species 
such  as  sea  turtles,  marlin  and  tuna,  seabirds,  and  cetaceans  (Davis  et  al.  2002,  Seki 
et  al.  2002,  Bakun  2006,  Polovina  et  al.  2006,  Yen  et  al.  2006).  Seki  et  al.  (2002) 
discussed  the  shift  in  the  area  of  highest  catch  in  the  1995  Hawaiian  International 
Billfish  Tournament.  During  this  particular  tournament,  the  area  of  highest  catch 
shifted  from  its  historic  location  to  the  region  coincident  with  the  convergence 
zone  at  the  edge  of  a  cyclonic  eddy  (Seki  et  al.  2002).  Additional  investigations  of 
eddies  in  the  lee  of  the  Hawaiian  Islands  show  them  to  be  areas  containing  greater 
concentrations  of  nutrients  and  phytoplankton  (Seki  etal.  2001,  Bidigare  etal.  2003, 
Vaillancourt  et  al.  2003).  Presumably,  this  leads  to  a  concentration  of  higher  trophic 
level  species  as  well. 

One  offshore  melon-headed  whale  track  was  associated  with  the  center  of  an  anti- 
cyclonic  eddy,  as  identified  by  the  significantly  slower  current  speeds  associated  with 
the  track  positions  (Track  3,  Fig.  2E— F).  Twelve  of  the  16  offshore  days  were  spent 
near  this  area  of  surface  convergence.  Several  studies  suggest  that  anticyclonic  eddies 
are  a  less  favorable  foraging  habitat  for  pelagic  marine  species  than  cyclonic  eddies 
(McGillicuddy  and  Robinson  1997,  Davis  et  al.  2002,  Bakun  2006).  However,  in 
their  study  of  seabird  associations  with  oceanic  eddies,  Yen  et  al.  (2006)  found  a 
correlation  between  several  species  and  both  cyclonic  and  anticyclonic  eddies.  The 
authors  attributed  this  correlation,  in  part,  to  the  aggregation  of  mid  trophic  level 
prey  such  as  fish  and  squid  that  occurs  in  the  convergence  zones  of  the  eddies  (Yen 
et  al.  2006).  Both  Howell  et  al.  (2010)  and  Polovina  et  al.  (2006)  reported  associ¬ 
ations  between  loggerhead  sea  turtles  (Caretta  caretta)  and  warm  core  anticyclonic 
eddies,  though  unlike  the  melon-headed  whale  track  discussed  above,  the  loggerhead 
positions  were  mostly  coincident  with  the  eddy  edges. 

Visual  inspection  of  the  five  offshore  tracks  also  identified  a  noneddy-associated 
area  southwest  of  Ni'ihau,  bounded  approximately  by  161. 5°W,  160°W,  20.5°N, 
and  22°N,  that  we  term  the  “Ni‘ihau  Box”  (Fig.  1).  No  clear  motivating  factor  could 
be  found  to  explain  why  these  animals  moved  through  this  area.  Seasonality  could 
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Table  2.  Details  on  melon-headed  whale  tracks  in  relation  to  the  Ni'ihau  Box  (see  text). 


Track 

Days 

offshore 

Days  inside  Ni'ihau  Box 
(%  of  offshore  time) 

Days  outside  Ni'ihau  Box 
(%  of  offshore  time) 

1 

11 

6(55) 

5  (45) 

2 

4 

2(50) 

2(50) 

3 

16 

2(12.5) 

14(87.5) 

not  be  adequately  investigated  as  the  three  tracks  spanned  only  4  mo  (April-July 
2008).  KS  tests  investigating  portions  of  the  tracks  inside  and  outside  the  Nihhau 
Box  separately  were  insignificant,  suggesting  that  neither  horizontal  nor  vertical 
currents  were  a  motivating  factor.  Bathymetry  was  also  ruled  out  as  all  offshore  track 
positions  were  in  waters  with  depths  greater  than  3,000  m.  However,  the  fact  that 
all  three  westward-moving  melon-headed  whales  independently  transited  this  area 
warrants  further  investigation.  Additionally,  of  the  three  melon-headed  whales  that 
transited  into  or  through  the  Ni‘ihau  Box,  two  spent  at  least  half  their  offshore  time 
in  this  area,  one  deviating  from  the  periphery  of  a  cyclonic  eddy  to  do  so  (Table  2, 
Track  1,  Fig.  2G-H). 

This  study  is  the  first  to  use  satellite  telemetry  data  to  gain  insight  into  the  pelagic 
habitat  of  melon-headed  whales.  Examining  melon-headed  whale  tracks  in  conjunc¬ 
tion  with  HYCOM  model  output  of  daily  SSH  and  sea  surface  currents  indicates  that 
when  melon-headed  whales  move  away  from  the  coastal  waters  of  Hawaih  they  move 
offshore  into  statistically  distinct  environments.  These  environments  coincide  with 
the  convergence  areas  of  both  cyclonic  and  anticyclonic  eddies,  indicating  that  they 
are  likely  selected  for  foraging  purposes.  Additionally,  all  westward-moving  melon¬ 
headed  whales  moved  into  or  through  a  small  area  southwest  of  Nihhau  for  reasons 
that  are  not  yet  apparent.  The  information  gained  through  this  study  may  also  indi¬ 
cate  pelagic  habitats  that  are  favorable  foraging  grounds  for  other  cetaceans  as  well. 
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